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FURTHER  STUDIES  ON  DYNAMIC  CRACK  BRANCHING 


M.  Ramulu,*  A.  S.  Kobayashl,*  B.  S. -J.  Kang,*  and  D.  B.  Barker** 


ABSTRACT 

The  newly  derived  dynamic  crack  branching  criterion  Is  verified  by  dynam¬ 
ic  photoelastic  analysis  of  dynamic  crack  branchings  In  thin  polycarbonate* 
single  edged  crack  tension  specimens.  Successful  crack  branching  was  observed 
In  four  specimens  and  unsuccessful  branchings  In  another.  Crack  branching 
consistently  occurred  when  the  necessary  condition  of  KIb  =  3.3  MPa  m»£nd  the 
sufficient  condition  of  r0  =  rc  =  0.7  mm  were  satisfied  simultaneously.  In 
the  unsuccessful  branching  test  the  necessary  condlton  was  not  satisfied  since 


Kj  was  always  less  than  K^. 


INTRODUCTION 


Crack  branching  represents  one  extreme  of  the  large  range  of  dynamic 
crack  propagation  behaviors  and  has  been  the  subject  of  numerous  theoretical 
and  experimental  Investigations*  several  of  which  can  be  found  In  References 
Cl* 23.  Recently*  the  authors  derived  a  crack  curving  and  a  branching  criteria 
based  on  the  directional  stability  of  a  propagating  crack  C3»4].  The  crack 
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curving  criteria  Is  a  micro-mechanical  model  of  continuous  micro-flaw  growth 
and  coalescence  In  the  vicinity  of  the  moving  crack  tip.  It  assumes  that  the 
crack  Is  momentarily  kinked  or  bifurcated*  when  an  off-axis  micro-flaw  con¬ 
nects  with  the  crack  tip  and  Is  the  dynamic  extension  of  the  crack  curving 

criterion  proposed  by  Strelt  and  Flnnle  CSoT 

//> 

The  dynamic  crack  curving  criterion  has  been  used  to  predict  the  crack 
kinking  angle  of  a  propagating  crack  under  pure  mode  I  as  well  as  mixed  mode 
condt1ons[33.  The  crack  branching  criterion  on  the  other  hand*  requires  a 
critical  stress  Intensity  factor  to  trigger  crack  branching  and  a  crack  curv¬ 
ing  criterion  for  predicting  the  crack  branching  angle  [4*63.  The  objective 
of  this  paper  Is  to  provide  further  evidence  In  support  of  the  dynamic  crack 
branching  criterion  advanced  by  the  authors. 


CRACK  BRANCHING  CRITERION 

The  crack  branching  criterion  [63  requires*  as  the  sufficiency  condition* 
a  crack  curving  criterion  [43.  The  latter  Is  based  on  the  postulate  that  the 
micro-cracks  ahead  of  the  crack  tip  dictate  the  direction  of  crack  propaga¬ 
tion.  When  an  off-axis*  l.e. *  9  t  0*  micro-void*  which  Is  within  a  critical 
distance*  r  ,  to  the  crack  tip*  is  actuated  by  a  critical  crack  tip  stress 
field*  It  deflects  the  crack  from  Its  otherewlse  self-similar  propagation 
path.  The  distance  between  the  micro-void  and  the  crack  tip*  rQ,  is  a 
characteristic  distance  which  Is  governed  by  the  singular  state  of  stress  as 
well  as  the  stress  acting  parallel  to  the  crack*  commonly  referred  to  as 
either  the  remote  stress  or  the  non-singular  stress  component*  oqx.  The 
critical  distance*  r  ,  Is  assumed  to  be  a  material  property. 
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The  angular  orientation  of  this  critical  micro-void  Is  determined  by  the 
maximum  circumferential  stress  criterion*  as  modified  by  Ramulu  and  Kobayashl 
C3L  which  assumes  that  the  crack  will  extend  towards  the  maximum  circumfer¬ 
ential  stress  at  a  distance*  r*  away  from  the  rapidly  propagating  crack  tip. 
Based  on  this  assumption*  under  pure  mode  I  loading*  l.e.,  K*  *  0.  Kn  =  0, 
the  condition  for  self-similar  propagation  of  a  straight  crack  Is  obtained  by 
setting  0  =  0  as 

(k. 


ro  =  128* 


where 


'I 


ox 


V0(C,C1*C2) 


Vo(c,c1c2)=B1  (c){-(l+s|)(2-3sf )-  ^1^(14+3S^)-16S1(S1S2)+16(1+S^)} 

1 


=  1  -  — 
1  2 

C1 


<•2  _  i 

S2  -  1  ‘ 

c2 


c  =  crack  velocity,  m/s 

c1  =  dilatational  wave  velocity,  m/s 

c2  =  distortional  wave  velocity,  m/s 


Here.  KT,  a  and  r  *  are  the  mode  I  stress  Intensity  factor*  remote  stress 

and  the  characteristic-  distance*  respectively*  and  can  be  determined  from  the 

current  dynamic  state  of  stress.  The  onset  of  crack  curving  of  a  rapidly 

propagating  crack  Is  governed  by  the  stability  of  the  propagating  straight 

crack  and  Is  assumed  to  occur  when  r  <  r  .  This  r  Is  a  characteristic 

o  c  c 

distance  derived  from  a  directional  stability  criterion  Involving  the  critical 
crack  tip  state  of  stress  where  0  suddenly  becomes  non-zero.  The  correspond¬ 
ing  angle.  0c*  for  a  maximum  circumferential  stress  can  be  determined  from  a 
transcendental  relation  Involving  the  critical  values  of  rc  and  ec  which  Is 
derived  from  maximizing  the  off-axis  maximum  circumferential  stress. 

•The  superscript  "dynw  to  Identify  dynamic  stress  Intensity  factor  will  not  be 
used  In  this  paper*  since  all  quantltltes  refer  to  dynamic  values. 
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The  crack  branching  criterion*  however*  Involves  not  only  the  critical  r£ 
but  also  the  maximum  Kj  as  a  necessary  condition  for  the  growth  of  multiple 
off-axis  secondary  cracks*  and  rQ  <  r£  as  a  sufficiency  condition  for  these 
multiple  cracks  to  kink  simultaneously.  Therefore*  the  crack  branching 
criterion  can  be  stated  as 

Kj  =  Max.  Kj  =  Necessary  condition 

r0  <  rc  Sufficiency  condition 

The  crack  curving  angle  0c  determined  from  the  latter  crack  curving  criterion 
Is  one  half  of  the  Included  crack  branching  angle. 

Crack  Branching  In  Homal lte-100  Fracture  Specimens  and  Pressurized  Pipes 

The  validity  of  the  above  crack  curving  and  crack  branching  criteria  was 
verified  by  dynamic  photoelasticity  results  of  Homal 1t-100  single  edged-notch 
(SEN)  specimens  and  the  wedge-loaded*  rectangular  double-cantilever  beam 
(WL-RDCB)  specimens  with  branched  cracks  [3*4*6].  Crack  branchlg  consistently 
occurred  when  the  dynamic  stress  Intensity  factor  reached  a  crack  branching 
stress  Intensity  factor  KIb  =  2.04  W*a/m  and  the  characteristic  distance  r0 
was  less  than  critical  distance  of  r  =  1.3  mm.  The  crack  branching  angles  of 
bifurcated  cracks  In  SEN  specimens  was  smaller  than  In  the  WL-RDCB  specimens. 
Differences  In  crack  branching  angles  are  expected  since  the  non-singular 
stress*  oox*  In  the  SEN  specimens  Is  compressive  and  suppresses  the  branching 
angle  whereas  the  tensile  0QX  in  WL-RDCB  enchances  the  branching  angle.  The 
crack  branching  angles  In  pressurized  steel  of  Reference  [7*8]  and  aluminum 
pipes  of  Reference  [9]  were  also  predicted  by  this  crack  branching  criterion 
[6]. 
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POLYCARBONATE  FRACTURE  SPECIMENS 

In  order  to  further  verify  the  dynamic  crack  branching  criterion#  a 
series  of  dynamic  photoelastic  fracture  experiments  Involving  thin  polycarbon¬ 
ate  fracture  experiments  were  conducted.  The  single  edged  notch  specimens 
with  blunt  starter  crack  were  127  x  225  mm  and  3.2  mm  or  6. 4  mm  thick.  At 
fracture  load#  the  crack  propagated  from  the  starter  crack  and  branched.  The 
dynamic  Isochromatics  surrounding  the  propagating  crack  were  recorded  with  a 
16  spark  gap  Cranz-Schardln  camera  system. 

The  fracture  parameters  of  KI#  Kjj  and  oQX  CIO]  associated  with  the 
running  crack  were  determined  by  least  square  fitting  a  theoretical  dynamic 
mixed-mode  stress  field  to  the  recorded  dynamic  Isochromatics.  The  Isochro- 
matlc  fringe  loops  were  digitized  and  analyzed  on  a  PDP-11/ 44  computer.  A 
least  square  algorithm  was  used  to  determine  KI'  KII  and  from  the 
multi-point  Isochromatlc  data  as  reported  in  Ref.  Cll#123.  The  estimated 
fracture  parameters  were  then  used  to  generate  the  corresponding  theoretical 
Isochromatics  which  were  superposed  on  the  experimental  Isochromatics  for  a 
visual  evaluation  of  the  accuracy  of  the  fitting  process.  A  flow  chart  of 
this  on-line  estimation  of  the  dynamic  fracture  parameters  from  the  recorded 
dynamic  Isochromatics  Is  shown  In  Figure  1. 


RESULTS 


Figure  2  shows  two  typical  dynamic  Isochromatlc  patterns  In  a  3.2  mm 
thick#  fracturing  polycarbonate  SEN  specimen.  At  fracture  load#  the  crack 
Initiated  and  propagated  from  a  blunt  starter  crack  with  several  unsuc- 
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cessful  attempts  at  branching  prior  to  a  successful  branching. 

Figure  3  shows  KI'  KII  and  <*ox  variations  associated  with  the  crack 
branching  experiment  of  Figure  2.  The  crack  after  Initiation*  propagated  with 
a  gradual  Increase  In  Its  dynamic  stress  Intensity  factor.  Immediately  prior 
to  branching*  the  Instantaneous  dynamic  stress  Intensity  factor  reached  Its 
maximum  value  of  3.3  MPa/m  with  negligible  KII  and  the  associated  remote 
stress*  aox*  attained  a  value  of  11.2  MPa.  By  smoothly  extrapolating  the 
average  Kj  and  KII  associated  with  the  two  branch  cracks*  an  after-branching 
Kj  *  2.2  MPavfii  and  KII  =  ±  0.  9  ►ff’ay'm  are  obtained. 

Figure  4  shows  two  frames  of  the  16-frame  dynamic  photoelastic  record  of 
a  propagating  and  branching  crack  In  a  3.2  mm  thick*  polycarbonate  SEN  speci¬ 
men.  The  crack  emanated  from  a  blunt  saw  cut  crack  and  propagated  through 
much  of  the  length  of  the  plate  with  Innumerable  unsuccessful  branching  prior 
to  the  successful  crack  branching.  Note  that  post-arrest  Isochromatics  sur¬ 
rounding  all  unsuccessful  branches  exhibit  a  pure  mode  II  crack  tip  deforma¬ 
tion. 


Figure  5  shows  the  dynamic  Kj,  Kjj  and  aQX  variations  obtained  from  the 
photoelastic  patterns  preceding  and  after  crack  branching  from  the  test  shown 
in  Figure  4.  Immediately  prior  to  the  crack  branching*  the  extrapolated  val¬ 
ues  of  Kj  and  0QX  at  the  onset  of  crack  branching  yielded  a  branching  stress 
Intensity  factor  of  Kj  =  3.3  MPa»4.  0QX  had  gradually  reached  a  value  of  11.5 
MPa*  which  Is  consistent  with  previous  test  results.  Immediately  after  bran¬ 
ching*  extrapolated  after-branching  the  average  Mode  I  and  Mode  II  stress  In¬ 
tensity  factors  of  Kj  «  2. 2  MPa>^m  and  Kjj  *  ±  0. 9  MPai4i  were  obtained. 
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Figure  6  shows  four  frames  out  of  a  16-frame  dynamic  photoelastic  record 
of  another  test  with  multiple  crack  branching.  Although  no  attempt  was  made 
to  analyze  these  post-branched  multiple  cracks*  the  data  up  to  the  onset  of 
successful  crack  branching  yielded  again  =  Kjb  =  3.32  MPavfiT  and  oqx  =  11.72 
MPa. 

During  the  last  ten  plus  years  of  dynamic  crack  branching  study*  we  have 

observed  that  either  the  unsuccessfully  branched  cracks  or  the  completely 

arrested  cracks  were  under  a  pure  mode  II  state.  In  the  tests  shown*  the  pure 

mode  II  Isochromatics  were  also  seen  (Figures  2*4*6)  at  the  unsuccessful 

branched  cracks.  Figure  7  shows  two  enlarged  views  of  Test  No.  KB-8208024  of 

Figure  5  with  a  mixed-mode  Isochromatlc  pattern  arres  of  the  branched  crack. 

Immediately  after  arrest  the  crack  tip  Isochromatics  transformed  Into  mode  II 

Isochromatics  with  the  45  y  second  Interval.  The  mode  I  and  mode  II  stress 

Intensity  factors  Kr  KII  and  oox  prior  to  and  after  the  crack  arrest  are:  Kj 

•  1.36  MPavfn,  KTT  =  0.06  MPavfii,  o  =  -7.7  MPa  and  KT  *  0.05  MPa»¥,  KTT  =0.7 
ll  ox  1  11 

MPa  m»  aQX  =  8.06  MPa  respectively.  This  suggests  that  the  arrested  branch 
crack  undergoes  a  mode  II  crack  tip  deformation  during  Its  unloading  loading 
process. 

Figure  8  shows  the  variations  of  KT„  «TT  and  o  associated  with  a 
straight  crack  with  unsuccessful  branches  In  3.2  mm  thick  specimen.  Although 
many  attempts  of  branchlgn  were  observed,  complete  branching  did  not  occur  In 
this  specimen  since  the  dynamic  stress  Intensity  factor  was  always  less  than 
the  Kjb  =  3.3  MPav'in.  Evaluations  of  two  additional  tests  yielded  the  crack 
branching  data  shown  In  Table  1.  The  critical  values  of  rc  ranged  from  0.6  to 
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The  crack  velocity  In  the  five  tests  were  essentially  constant  and  at 
about  23  percent  of  the  dllatatlonal  wave  velocity*  c^  =  1955  mps.  The  same 
velocity  was  observed  In  the  crack  curving  experiments  conducted  with  this 
material  C13L  It  appears  that  this  crack  velocity  Is  the  maximum  observed  In 
all  the  dynamic  fracture  tests  Involving  polycarbonate. 

Figure  9  shows  the  variations  In  characteristic  distance  rQ,  which  was 
computed  by  Equation  (1)*  for  the  propagating  cracks  prior  to  the  crack 
branching  In  the  five  tests.  Although  the  value  of  rQ  has  a  scatterband  of 
0.5  to  0.9  mm*  as  shown  In  Figure  8*  all  extrapolated  rQ  at  crack  branching 
reached  an  average  minimum  value  of  0.7  mm.  This  crack  branching  rQ  =  rc  = 
0.7  mm  represents  the  sufficiency  condition  for  crack  curving  and  Is 
consistent  with  the  rc  value  determined  from  the  crack  curving  experiments  of 
polycarbonate  material  [133. 

Table  1  shows  the  crack  branching  stress  Intensity  factors*  KIb,  the 
critical  distance*  rc,  and  the  measured  and  predicted  crack  branching  angles 
of  all  four  test  results  of  successful  crack  branchings.  The  dynamic  stress 
Intensity  factor  at  the  onset  of  crack  branching  reached  an  average  maximum 
value  of  3.3  MPa*^n.  This  branching  stress  Intensity  factor  was  found  to  be 
independent  of  the  thickness  of  the  specimen  as  well  as  the  Initial  and 
branching  crack  lengths. 


Crack  branching  angles  were  computed  by  using  the  crack  curving  criterion 
and  are  listed  In  Table  1.  The  average  crack  branching  angle  Is  25°.  The 
branching  angle  In  this  series  of  SEN  specimens  varied  between  22°  and  34°  and 
Is  consistent  with  our  previous  results  Involving  Homallte-100  [4]. 

DISCUSSION 

Post-branching  cracks  In  all  tests  always  curved.  Kalthoff  [14], 
observed  that  the  direction  of  two  branched  cracks  attraction  or  repulsion.  Is 
controlled  by  The  photoelastic  patterns  of  running  branched  cracks 

showed  that  the  crack  was  perpendicular  to  the  load  direction.  This  strongly 
suggests  that  the  crack  runs  parallel  to  the  compressive  stress  direction  even 
under  mixed  mode  conditions  which  exist  after  branching.  Therefore  post¬ 
branching  crack  propagation  Is  also  strongly  dependent  on  the  Kjj/Kj  ratio  as 
well  as  on  a  .  Figure  10  shows  the  post-branching  crack  curving  of  specimen 
No.  820822.  The  measured  and  calculated  angles  are  marked  on  Figure  10  show 
that  the  crack  curving  angle  gradually  decreased  In  magnitude  along  with 
Increase  In  negative  aQx  and  Is  In  agreement  with  the  numerical  results  of 
Ref.  [153. 

Figure  11  shows  the  typical  fracture  surface  In  a  6.  4  mm  thick  polycar¬ 
bonate  SEN  specimen  associated  with  the  crack  branching.  Clear  mirror,  mist 

and  hackle  zones  are  visible.  This  fractured  surface  Indicates  that  the  o 

ox 

term  which  Is  the  paral lei  stress  in  this  specimen.  Is  compressive  and  opens 
the  micro-cracks  In  the  form  of  tongues.  Although  the  dynamic  stress 
Intensolty  factor  Is  almost  equal  to  Kjb,  the  crack  did  not  branch  at  the  fine 
hackle  zone  but  branched  when  this  hackle  zone  became  rougher  and  when  the 
sufficiency  condition  was  met. 


CONCLUSIONS 


1.  Dynamic  crack  branching  criterion  proposed  by  the  authors  successfully 
predicted  the  crack  branching  when  the  necessary  condition  of  K,  *  K.. 
which  triggered  the  generation  of  secondary  cracks*  and  the  sufficiency 


condition  of 
satisfied. 


<  V 


which  caused  the  crack  tip  diversion*  were 


2.  A  crack  branching  stress  Intensity  factor  of  K-.  *  3.3  MPavfo  and 

characteristic  radius  of  r  *  0.7mm  are  determined  for  this  polycarbonate 
sheet. 


3.  Crack  curving  of  post-branched  cracks*  attraction  and  repulsion,  depends 
not  only  on  Kj./K,  but  more  Importantly  oncf  .  Negative  a  suppresses 
the  crack  curving1 Irrespective  of  the  sign  9  Kn/K r 
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